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Abstract 

We derive the large order behavior of the perturbative expansion for the continuous 
model of tethered self-avoiding membranes. It is controlled by a classical configu- 
ration for an effective potential in bulk space, which is the analog of the Lipatov 
instanton, solution of a highly non-local equation. The n-th order is shown to 
have factorial growth as (— cst) n (n!)( 1-e/ '' D ), where D is the "internal" dimension 
of the membrane and e the engineering dimension of the coupling constant for self- 
avoidance. The instanton is calculated within a variational approximation, which is 
shown to become exact in the limit of large dimension d of bulk space. This is the 
starting point of a systematic 1/d expansion. As a consequence, the e-expansion 
of self-avoiding membranes has a factorial growth, like the e-expansion of polymers 
and standard critical phenomena, suggesting Borel summability. Consequences for 
the applicability of the 2-loop calculations are examined. 
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1 Introduction 



Flexible polymerized tethered membranes (also called polymerized membranes or tethered 
membranes) exhibit fascinating statistical properties [1,2]. Tethered membranes with only 
short-range repulsive interactions ("self-avoidance") may be seen as the 2-dimensional 
analog of polymer chains in a good solvent. They are expected to exist in two phases: 

(1) a flat (or rigid) phase with an infinite persistence length but non-classical elastic 
properties, not found for polymers which have a finite persistence length [3-5], and (2) a 
crumpled phase similar to that of polymers [6,7]. At variance with polymers, the crumpled 
phase is always swollen (whatever the dimension d of space is) and has a configuration 
exponent v larger than the mean-field exponent z/ mf = 0, whereas for polymers it is well 
known that for d > 4 one has v = z/ mf = 1/2 while for d < 4 one has v > z/ mf . This follows 
from a simple dimensional argument, which shows that self- avoidance is always relevant 
(in the renormalization group sense) [6]. 

The properties of the crumpled swollen phase have been much studied theoretically. 
In particular, a systematic framework for renormalization group (RG) calculations has 
been proposed in [8-10]. It consists in: (1) a self-avoiding membrane model (SAM) which 
is a non-trivial extension of the continuous Edwards model for polymers [11,12], and 

(2) a renormalization procedure which generalizes the direct renormalization [13] for the 
Edwards model. For 2-dimensional tethered membranes, the upper critical dimension d uc 
is infinite, but the SAM model can be extended to membranes with non-integer internal 
dimension D < 2, which have a finite d uc . Then the framework of [8,9] allows to perform 
a double expansion both in D (the internal dimension of the membrane) and d (the 
dimension of bulk space) to study the physical case (D = 2, d < oo) of a 2-dimensional 
membrane in finite dimension d. The starting point for the e-expansion (characterized 
by the dimension Do < 2 for the membrane) can be chosen arbitrarily. This provides an 
additional parameter for the expansion, which can be used to optimize the calculations 
[14,21]. 

This approach, which amounts to perform perturbative renormalization group calcu- 
lations for .D-dimensional membranes, has raised considerable challenges. Whereas for 
polymers (i.e. when D=l) the Edwards model can be mapped onto the local $ 4 field the- 
ory with n = components [15], and direct renormalization is equivalent to the standard 
renormalization (in the MS scheme) [16], this equivalence does not hold for D ^ 1. It 
was finally shown in [17,18] that the perturbative RG calculations are mathematically 
consistent, by proving that the SAM model is renormalizable to all orders in perturbation 
theory. If perturbative calculations are simple at first order (with some subtle points [19]), 
they present considerable difficulties at second order, and require a lot of analytical and 
numerical work [20,21]. 

An important issue is to understand if these calculations make sense beyond perturba- 
tion theory, or if non-perturbative effects destroy the consistency of the approach. A first 
step is to understand the large order behavior of perturbation theory for the SAM model. 
To our knowledge, nothing was known about this problem up to now, except for D = 1, 
where one can use the equivalence to the n — > limit of <£> 4 -theory, for which a non-trivial 
solution of the equation of motion at negative coupling constant, the so-called instanton, 
governs the large order behavior [22,23]; however, this analogy does not provide a physical 
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picture of which polymer-configurations dominate the large orders. In this paper we shall 
formulate the problem of the large order behavior for the Edwards model, in a way which 
is directly applicable both to polymers and to membranes. 

Let us summarize the organization and the results of this paper. In Sect. 2 we recall 
the definition of the SAM model. Then we develop a general semiclassical argument 
to compute the large orders of the perturbative series for the SAM model. Using the 
formulation of the SAM model as a model of a "phantom" membrane (without self- 
avoidance) in a random imaginary external potential V, we show that the large orders 
are controlled by a real classical configuration for this potential V, which is the analog 
for SAM of the instanton for <3> 4 . This "SAM instanton" potential V is the extremum of 
a non-local functional «S[V], which cannot be calculated exactly. We obtain the general 
form for the asymptotics of the term of order n, which is 

n d/* (_ C )» ( n |)WJ> , (1) 

where D is the internal dimension of the membrane, d the dimension of bulk space, 
e = 2D — d(2 — D)/2 the engineering dimension for the self-avoidance coupling in the 
SAM model and C a positive constant depending on D and e (or d). This behavior is 
universal: the constant C obtained from the instanton does not depend on the internal 
shape or topology of the membrane. 

In Sect. 3 we show that for polymers (D = 1) our results fully agree with the classical 
results on the large orders for the $ 4 theory, as obtained from instanton calculus. We 
then give a simple physical interpretation of our instanton for the SAM model: the SAM 
instanton describes the metastable equilibrium configuration of a membrane submitted 
to two competing forces, a contact attractive interaction and a global entropic repulsion 
due to thermal fluctuations. 

The SAM instanton equations are in general not solvable. In Sect. 4 we propose a 
Gaussian variational approximation scheme and find the instanton potential V and the 
large order constant C within this approximation. By construction, this is a lower bound 
on the exact result. 

The variational results are discussed in Sect. 5. They qualitatively agree with the 
exact results for polymers (D — 1). We estimate the "optimal order" for the e-expansion, 
and we argue that the estimates for the configuration exponent v obtained from second 
order calculations [20,21] are reasonable, especially for large bulk dimension d. Finally we 
observe that they imply that the perturbative series for the SAM model should become 
"quasi-convergent" in the limit d — > oo. 

In Sect. 6 we study more carefully the validity of the variational approach. We show 
that it becomes exact in the limit of large space dimension d — > oo. Technically one has 
to keep e fixed, that is to let the internal dimension D of the membrane go to 2. To 
go beyond the leading order, we construct a systematic perturbative expansion around 
the variational solution, and show that this expansion can be organized as a systematic 
expansion in powers of 1/d, by resumming infinite classes of diagrams. We compute 
explicitly the first correction in 1/d for the instanton V and the large order constant C. 
In particular, we find that along the critical line e = 0, the instanton potential V has 
log(e) singularities, but that the large order constant C is finite. 

In Sect. 7 we present our conclusions and discuss open problems. 
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2 Large orders and instantons for the SAM model 



2.1 The SAM model 

We consider a D-dimensional manifold M. with size L and volume V = LP (typically the 
D-dimensional torus T D = [0,L] D ) in <i-dimensional Euclidean bulk space. A configura- 
tion of the manifold is described by the continuous function 

x eM ^ r(x) e R d . (2) 

The partition function is 



Z 

with the Hamiltonian 



[b;L] = Jv[r\ e - n ^ (3) 



W\ b, L] = J d D x l -{Vr{x)f + b -J L d D x J d D y 5 d (f(x) - f(y)) . (4) 

b > is the repulsive 2-point interaction coupling which describes self-avoidance. The 
functional integration measure X>[f] = Yi x <i d r(x)/Z Q is normalized such that the partition 
function of the free Gaussian manifold Z[b — 0, L\ — 1. The canonical dimensions of x, f 
and b are 

H = -i , H = ^ . H = 2D -^ = £ . ,5) 

We know that the partition function Z and expectation values of in bulk space trans- 
lationally invariant operators (O) are well defined as a perturbative series in the coupling 
constant b, as long as the dimension of the manifold D is less than 2 and as long as e is 
positive 

< D < 2 , e > (6) 

and provided that the size V of the manifold is finite 1 . For e > 0, perturbation theory 
suffers from strong Infra-Red (IR) divergences when the size of the manifold becomes large 
L — > oo, which signal a breakdown of mean-field theory and the appearance of anomalous 
dimensions in the scaling properties of large self-avoiding manifolds. 

For e — > 0, physical UV singularities appear at all orders in perturbation expansion. 
They can be absorbed in a renormalization of the coupling constant b and a rescaling of 
bulk space (wave function renormalization) 

b = b R Z b (b R ) , r = Z l l\b R )r R (7) 

so that physical observables are UV finite at e = when expressed in terms of renormal- 
ized quantities b R and r R . The theory, although non-local in the internal a>space of the 



1 Strictly speaking, for < e < D, the first terms in the expansion of the partition function suffer from 
short-distance divergences. These Ultra- Violet (UV) singularities are finite in number and can be recast 
in a counterterm proportional to the volume V of the manifold. Moreover they disappear in physical 
observables; they are thus of no physical significance and unimportant for the purpose of this paper. 
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membrane, is renormalizable [17,18]. As a consequence, it is possible to write renormal- 
ization group equations which encode how the effective renormalized coupling flows with 
the length scale, and in particular with the size L of the membrane. The large L behavior 
of the membrane is governed by an IR stable fixed point b* R , which can in principle be 
calculated order by order in perturbation theory as a series in e, similar to the celebrated 
Wilson-Fisher e-expansion. In practice, such calculations are very difficult, and have been 
performed only up to second order [20,21]. 



2.2 Instanton calculus 



By dimensional analysis, the partition function (3) only depends on the dimensionless 
coupling constant 

9 = bU (8) 



via 



and is defined as a series 



Z[b;L] = Z[g,L=l] = Z[g] 



Z[g] = X>«S n 

n=0 



(9) 
(10) 



Of course, Z[g\ also depends on the shape of the manifold Ai. 

Let us assume that Z[g] is analytic around the origin for — n < arg(g) < ir, and has 
a discontinuity along the negative real axis. This assumption is natural, since for g < 0, 
the membrane is in a collapsed state and the perturbation expansion is performed around 
an unstable classical state. Then we can write dispersion integral 



— g- n - l Z[g] 



f * g—Hm(Z[g + *0 + ]) . (11) 
Jo 71 



To obtain the behavior for large n, it turns out that it is sufficient to evaluate the integral in 
(11) in a saddle point approximation. Indeed, we shall show that, at least for < e < D, 
the integral is at large n dominated by the discontinuity of Z[g] at small negative g. 
Moreover, Z[g] is dominated by a saddle point, when re-expressed as a functional integral 
over properly defined auxiliary fields. 

The Hamiltonian (4) is non-local and involves a distribution of the field f. It is 
convenient to introduce as an auxiliary field the density of the membrane in bulk space 

p(f) = J d D x 5 d (f-f(x)) 

and to write the interaction term as 

d D y5 d (f(x) -r(y)) = 



(12) 



d D x 



dVp(f)' 



(13) 



Introducing a potential V{r) conjugate to p(r), we can re-express the interaction term in 
the partition function through a Hubbard-Stratonovich transformation as 



exp 



- I d D x 



d D y5 d (f(x)-f(y)) 
V(r[ 



J V\p\ J V[V] exp J d d r \ V(r) p{r) - J d D x 5 d {f- f{x)) 



(14) 
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where 

/ v ^ = f +0 ° n Mr-) , / v[v\ = f 100 n . ^ 

J J-oo r - GRd J J -ice ^ 

Inserting Eq. (14) into Eq. (3) and integrating over p yields 

Z(b;L) = Jv[r\ Jv[V]e- H ^ v ' b ^ (16) 

with the new effective Hamiltonian 

n'[f,V;b,L] = J d D x Q(Vf(x)) 2 + V(r%x)) y j - ^ J d d f V(f) 2 . (17) 

This representation is nothing but the generalization of the well known formulation of 
the Edwards model as a model of free random walks in an (imaginary) annealed random 
potential. As above, Z is a function of the dimensionless coupling g and we replace b — > g 
and L — > 1 as in Eqs. (8) and (9). 

As argued before, we aim in calculating the partition function for small negative g. 
For that purpose, it is convenient to rescale the coordinates and the potential V(r) 

x^(-g)T^x , f^(-g)^f , V^(-g)^V, (18) 

so that we now consider a membrane with size L = {—g)^~ € and volume 

V = L D = (-g)T^ . (19) 

This yields the rescaled Hamiltonian 

K*c[r,V;L] = f_d D x Q(Vr») 2 + V(r»)) - V - J d d fV(f) 2 . (20) 

The integral over f for fixed potential V defines the free energy density £[V] of a "phan- 
tom" (i.e. non self-avoiding) membrane in the external potential V 

Z[V] = e- 9£ [ y l = Jv[r\e-k dDx ^ x))2 + v W x)) . (21) 

The partition function finally becomes 

Z(g) = JvlV]e-^lv]+y^rV(ry} (22) 

The crucial point of this formulation is that according to Eq. (19), as long as 

< e < D , (23) 

the limit g — > _ corresponds to the thermodynamical limit when the volume V — > oo. 
In this limit the free energy density £[V] has a finite limit (from extensivity) so that the 
volume appears only as a global prefactor in the exponential of (22). Hence in the large 
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V limit the integral (22) is dominated by a saddle point V^ nst , which is an extremum of 
the effective energy S[V] for an infinite and flat membrane. The latter is defined as 



S[V] = S[V\ + \J d d fV(r) 2 , 



(24) 



where S[V] is defined in (21) as the free energy density of an infinite flat membrane in 
the potential V. This saddle point Vi ns t(r) is the non-trivial instanton, since the action S 
of the trivial extremum V(r) = is real and does not contribute to the discontinuity of 
2(g). Moreover, as the instanton is obtained through the thermodynamical limit L — > oo, 
it is independent of the shape of the initial membrane. This implies that the large order 
behavior of perturbation theory is universal, and does not depend on the internal geometry 
of the membrane. Let us now derive the saddle-point equations: 
The variation of the free energy density is in general 



5£[V] 
SV(r) 



(5[r\) 



v 



(25) 



where S[r\ is the normalized density of the membrane 



5[r\ 



v 



y J d D xS d (f-f(x)) 



(26) 



(which has a finite limit when the volume becomes infinite), and ( ) v denotes the expec- 
tation value for the phantom membrane in the potential V, as defined in Eq. (21). Hence 
extremizing S[V] leads to the variational equation for the instanton potential V inst 



= {S[r])v inBt + V inst (r) 



(27) 



Let us postpone the solution of Eq. (27) and first ask what the consequences of the 
existence of an instanton for the large order behavior are. Denoting by <Si nst the action 
for the instanton <S[Vj nst ], we deduce from Eqs. (19) and (21) that for small negative g, 
the discontinuity of 2(g) behaves as 



lm(2(g)) « exp S 



'inst 



(28) 



and the integral representation for z n (11) can be evaluated by the saddle point method 
at large n. This saddle point is at 



9c = 



c. 



inst 



n(l-e/D) 



l-e/D 



(29) 



and replacing the integral in (11) by its value at g c gives the large n behavior at leading 
order 

l-e/D 



z n ~ (-C) n (n\)^/ D 



l-e/D 
Si 



inst 



(30) 



Let us briefly discuss this result. For < e < D, perturbation theory is divergent with 
alternating signs. For e = 0, one recovers the typical factorial behavior (— C) n n\ of field 
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theories, provided that S mst remains UV finite. As we shall see in the next section, our 
result (30) coincides for D = 1 with the large order behavior deduced from the $ 4 model 
with n = components. The reasoning seems to break down at e = D, but we shall see 
that in fact the factor of C, when considered as a function of D and e, is regular at e = D 
and can be continued to the region e > D. Thus, the asymptotics (30), although derived 
for < e < D, is valid in the whole physical domain < e < 2D. A more rigorous 
argument is as follows: Eqs. (28) and (29) are still valid for e > D; the instanton then 
governs the behavior of the discontinuity of Z(g) at large g. This means that the saddle 
point of Eq. (30) for large n now is at large negative g. 

To go beyond this estimates, one must (i) compute the instanton and its action, and 
(ii) integrate the fluctuations around the instanton in (22). If one assumes that this 
calculation goes along the same lines as in standard field theory, one must first isolate 
the zero modes, i.e. the collective coordinates of the instanton. As we shall see later, the 
instanton Vi nst is rotationally invariant and is characterized by its position in rf-dimensional 
space only. Thus it has d zero modes, each of them gives a factor V 1//2 (by a standard 
collective coordinates argument), and the remaining fluctuations 5±V (orthogonal to the 
translational variations S^V ~ * ) give a finite determinant *A.. Therefore we expect 
the semiclassical estimate for the discontinuity to be 

lm(Z(g)) ~ ^-V2 y| e -V5 inst ( 31 ) 
and that the large n behavior is more precisely 

z n ~ A'n d/2 {-C) n (n!) 1 ^ [1 + . . .] . (32) 

Finally we shall see that the action of the instanton remains finite in the limit e — > 0. As 
in standard $ 4 theory, one expects UV divergences to appear only for fluctuations around 
the instanton, and that these divergences are cancelled by the same renormalizations as 
in perturbation theory. This would imply that our large order estimate (30) is also valid 
for the renormalized theory at e = 0, in particular for the renormalization group functions 
which enter into the e-expansion of the scaling exponents. Renormalization however has 
to be taken into account when evaluating the constant A! in (32). 



3 The polymer case and physical interpretation of the instanton 

Before discussing membranes, let us study in detail the special case D — 1, where the 
model reduces to the Edwards model for polymers. Using the well known mapping be- 
tween the problem of a Brownian walk in a potential V(r) and quantum mechanics of a 
single particle in the same potential, the free energy density £[V] of a linear chain fluc- 
tuating in a potential V(r) is in the thermodynamic limit given by the lowest eigenvalue 
E of the operator 

H = -j + V(f) , (33) 
where A is the Laplacian in d dimensions. Thus we have 

£[V\ = Eq . (34) 
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Denoting by ^o(f) the ground state wave-function, and using Eq. (27) and the standard 
result from first order perturbation theory 

we obtain the instanton potential 

V inst (f) = - (^o(rl) 2 . (36) 
The eigenvalue equation H^/q = Eq^q becomes non-linear 

^Af + £„f + V = . (37) 
Since obeys the normalization condition 

ll^oll 2 = J d d f^ (r) 2 = 1 , (38) 

the wave function ty Q and the ground state energy E are fully determined by Eqs. (37) 
and (38). Eq. (37) has nontrivial normalizable solutions for 2 < d < 4 and Eq < 0. In 
addition, the ground state ^q is rotational symmetric, i.e. does not vanish at finite f. The 
action for the instanton (24) finally reads 

<S inst = E + l - J d d f * 4 . (39) 

To make contact to the instanton analysis in the Landau-Ginsburg- Wilson (LGW) 
$ 4 -theory with n = components, remark that Eqs. (37)- (38) hold if and only if ^ and 
E are extrema of the action 



S'[V,E} = E + J d d f 



1 1 
_(V^) 2 - E^ 2 # 4 

2 y 1 ' 2 



(40) 



This is the standard Landau Ginsburg Wilson action with negative coupling associated 
to \l/ 4 and mass m 2 = —IE. Moreover, at the extrema, the two actions are equal 

<S inst [*o,£o] = S'[* ,Eo] . (41) 
The relation becomes clearer by the change of variables 

*« = (izf) 1 ■ < 42 > 



The action S' then reads 



-2E\ 2 ^ 

S'[V,E] = E + [—A 5 LGW [$] (43) 



9 




with 



>LGW 



m = / 



I( V $) 2 + L^$2 _ I$4 

2 K ' 2 2 



(44) 



We can extremize (43) with respect to E and <3> independently, and denoting by <3>o and 
Eq these extremizing solutions, we get 

E = d-2)S LGW [%]^ . (45) 

The change of variables in Eq. (42) was constructed such that the instanton action takes 
the simple form 

S iQSt = S'[^ ,E ] = {I - 1) S LGW [%]^ . (46) 

Since for polymers (D = 1) d/2 — 1 = 1 — e/D, we can use Eq. (30) to write the large 
order constant C of the Edwards model as 

i = S LGW [%] . (47) 

This result could have been derived directly from the standard field theoretical formulation 

of the Edwards model as a n = component ($ 2 ) model. 

The equation for the instanton derived from the action (44) admits a regular solution 
$o(M) f° r an Y < d < 4, so that nothing special occurs at the point d = 2 (i.e. e — D — 1) 
as one might have expected from Eq. (30). Let us note that since the "mass" in Eq. (44) 
is equal to 4 — d, it is positive for d < 4 but vanishes at the critical dimension d — 4, so 
that the instanton solution $ still exists for d = 4. In Fig. 1 we plot C^ 1 (d) for < d < 4, 
as obtained from numerical integration. 

It is interesting to give a physical interpretation of the instanton for the Edwards 
model, since this interpretation is the same for membranes with D ^ 1. Let us first recall 
the standard interpretation of the instanton for the LGW model with action (44), i.e. 
negative $ 4 -coupling. The classical false vacuum \I/(r) = is separated from the true 
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vacua \J/(r) = =Foo by a finite barrier. The instanton solution ty Q describes a metastable 
droplet of true vacuum (with \l/o( r ) 7^ inside the droplet) in the false vacuum, which is on 
the verge to nucleate. Indeed, if the droplet is slightly larger, the positive surface energy 
dominates and the droplet shrinks and finally vanishes, while if it is slightly smaller, the 
negative volume energy dominates and the droplet expands. 

We now consider the energy density S[V] given by Eq. (24). It corresponds to the 
total free energy of a polymer globule trapped in the potential well V(r) < 0, where 
this effective potential results from the attractive 2-point interaction between elements of 
the polymer (since we are at negative coupling, b < 0). To see how S varies with the 
average gyration radius of the polymer, it is convenient to consider the following scale 
transformation on V 

2D 

V(r) -> V x {r) = A^ V(Xr) . (48) 
Simple dimensional analysis shows that under (48) 

S[V] -> \^£[V] , Jv 2 ^\^Jv 2 (49) 

(here D = 1 and e = 2 — d/2). As long as e < D, and for large A, i.e. when shrinking 
the polymer globule, it is the first term £ < which dominates and the total free energy 
S becomes large and negative; while for small A, i.e. when expanding the globule, it is 
the second term on the r.h.s. of (24) which is larger than 0, and which dominates. Thus 
in this mean field picture, i.e. neglecting thermal fluctuations around the instanton, large 
globules tend to expand, while small globules tend to collapse. This has a simple physical 
interpretation: the polymer trapped in its own potential is subject to two opposite forces, 

(i) attractive forces between its elements which would like to make the polymer collapse, 

(ii) entropic repulsion which exerts a pressure on the well and would like to expand the 
polymer (until it becomes a free random walk). What our calculation implies is the simple 
fact that for large radius (i.e. small A) entropic repulsion dominates, while at small radius 
(large A) attraction dominates and the polymer collapses. Thus the instanton solution 
describes a polymer with attractive interactions on the verge to collapse into its dense (and 
most stable) phase; this is similar to the instanton in the LGW theory which describes a 
bubble of true vacuum on the verge to nucleate and to destroy the false vacuum. 



4 Gaussian variational calculation 

For D^l (and in general for < D < 2 non-integer) we know of no exact method to 
calculate the instanton. A simple and natural approximation is the variational method, 
i.e. the Hartree-Fock approximation. 

To evaluate the free energy density £[V] of the free, i.e. non-interacting membrane in 
a potential V, and described by the Hamiltonian 




we introduce the trial Gaussian Hamiltonian 

n var = I d D x I d D y±r(x)K(x-y)r(y) 
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d D k 1 



f(k)K(k)f(-k) , 



(2n) D 2 

where ~ denotes the Fourier transform. The free energy for the trial Hamiltonian is 



(51) 



£var = - y Log 



V[r\ e 



(52) 



and the factor of 1/k 2 comes from the normalization of the measure T>\f\ taken such that 
£[V — 0] = 0. V is the total volume of the membrane. 

The Hartree Fock approximation amounts in replacing £[V] by the best variational 
estimate £ var [V] 

£[V] < £ vat [V] = £ vaT + i (H v -H yar ) vax . (53) 

( )var denotes the average with respect to the trial Hamiltonian 7i var and one must look 
for the trial Hamiltonian 7i var (i.e. the kernel K) which minimizes £ V ar|X]- Denote by 
V(p) the Fourier transform of the potential V(r). Since the variational Hamiltonian is 
Gaussian, it is easy to compute the second term on the r.h.s. of (53), V~ 1 {Hy — 7^var)var 
in the infinite volume limit: 



(V(r(0)) 

d d p 



var + 2 



(Vf(0)) 2 ), 



d D xK(x) (f(x)f(O)), 



V(p) (e 



ipr(0)\ 



(2vr) d 



+ 



d D k e -R(k) * * 

— (r(k)r(-k)}, 



-Q 



(2tt) d 
d D k 1 



2 J (27r) D K(k) 



d D k 



d 

' 2.1 {2tc) d \K(k) 



(- 



k' 2 



(54) 



Combining Eqs. (24), (52) and (54), we finally obtain the variational estimate for the total 
energy of the instanton 



S var [V] = £ var [V] + l -J d d fV( 
d d p 



(27T) 



d 

+ 2 



d D k 
(2^p 



Log 





r d D k i \ 




W D k(k)J 


K{k) 
k 2 


+ e -i) 

K(k) J 



+ \v{p)V{-p) 



(55) 



We now extremize Eq. (55) both with respect to K (variational approximation) and with 
respect to V (to obtain the instanton solution). Extremizing w.r.t. K{k) yields the equa- 
tion 



k(k) 



k - 



d d p 



d J (2ir] 



d 0*V(p) exp 



p 2 



d D k 



2 J (2n) D K(k) 



(56) 
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which implies that the variational Hamiltonian depends just on a mass m v 

K(k) = k 2 + ml r . 
Extremizing Eq. (55) w.r.t. V{p) gives 



with 



A 



d D k 



exp 



■— A 



m 



D-2 
var 



r(i-f) 



(57) 



(58) 



(59) 



{2n) D K{k) "" var (4tt)^/2 

T is Euler's Gamma function. Thus, in the variational approximation the instanton 
potential is Gaussian. Inserting Eq. (58) into Eq. (56) yields the self-consistent equation 
for m var 

d d p 



(2tt)^ ~ 2 
We finally get in terms of D, e and d = 2{2D - e)/(2 - D) 



2 r(--)i+s 



d\ D- 



The final result for A reads 



1 „ /Q nN -' 2 D-2 
' "Z-D\ D-e 2^37 



(60) 



(61) 



(62) 



We can now insert these results into Eq. (55), and after straightforward calculations get 
the variational instanton action 



(63) 



The corresponding variational estimate for the large order constant C defined by Eq. (30) 
is 



(64) 



As claimed in the previous section, although intermediate results are singular at e = D, 
the final result is regular for all e > 0. We shall discuss the physical significance of these 
results in the next section. 



5 Discussion of the variational result 

5.1 D = l 

It is interesting to compare the variational estimate with the exact result for polymers, 
i.e. for the case D — 1. Let us consider the LGW instanton action, as given by Eq. (44). 
It is equal to the inverse of the large order constant C. On Fig. 1 we have plotted the 
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variational result for 1/C var , as given by Eq. (64) and the exact result for 1/C obtained by 
numerical solution, as a function of < d < 4. First we note that always 

C > C var (65) 

as expected from the variational inequality £ < £ v3iT . This implies that the variational 
method gives an underestimate of the large orders. 

Second the variational estimate becomes good for small d, and exact for d — > 0. This 
is not unexpected, since in that limit the membrane Ai has no inner degrees of freedom, 
and the functional integration over V(r) reduces to a simple integration over V e M. 
Since this integral is Gaussian, the variational method becomes exact. 

Finally, the variational estimate for C is regular when d — > 4, and then equals l/(2n 2 ); 
this is 50% smaller than the exact result 3/(47r 2 ). Thus the variational method is only 
qualitatively correct when e = 0. This is not so surprising, since the limit e — > is 
somewhat peculiar. Indeed when d — 4 the ground state energy E Q in the equation (37) 
for the wave function is then equal to 0. Then the most general solution to Eq. (37) 
(for d — 4 and E = 0) is 

*o(r) = -r-^ , (66) 
7q + r 2 

with tq an arbitrary scale (the size of the instanton). tq is fixed by the normalization 
condition (38) which cannot be fulfilled at d — 4 for finite tq. In fact a more careful 
analysis of the rotationally invariant solutions of Eqs. (37) and (38) (see Appendix A) 
shows that as d — > 4, _E should scale as _E ~ 4 — d and that for < 4 — d <C 1 the true 
solution is well approximated by Eq. (66) (at least as long as |r| 2 (4 — d) <C 1) with an 
instanton size r which vanishes as d — > 4 as 

r ~ — = . (67) 

V|log(4-d)| 

The corresponding instanton potential Vi nst = — l^ol 2 is also singular in the d — > 4 limit 
(it may be considered as a Dirac-like 5-function), and is very poorly approximated by the 
Gaussian variational solution at D — 1, d — 4 for the potential 

\/ var = - 16vr 2 e- 16 " 2 " 2 , (68) 

with positive width. As usual with variational methods, the approximation for the ground 
state energy is much better than that for the wave function. 



5.2 Consequences for the e-expansion 

Of course, one is interested in the consequences of these large order estimates for the 
e-expansion of the scaling exponents for self-avoiding membranes and polymers. Let us 
recall that in renormalized perturbation theory one computes the renormalization group 
/3-function /3(g), as a power series in g of the form 2 

(3(g) - - eg + B l9 2 + 0(g 3 ) . (69) 

2 Strictly speaking g is now the renormalized coupling constant 
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Its zero at g* = e/B 1 + 0(e 2 ) is the IR fixed point which governs the scaling limit for 
large membranes. Other anomalous dimensions, like the dimension u(g) of the field f 
(which gives the fractal dimension of the membrane) can also be computed as a series in 
g. Their values at the fixed point g* give the scaling exponents of the membrane, and 
may be expanded as power series in e. 

By analogy with the ordinary Wilson-Fisher e-expansion for LGW field theories, let us 
assume that the large orders of the function (3(g) and of the other anomalous dimensions 
are given by the instanton estimate, and that they can be resummed by Borel-techniques. 
We are not able at the moment to give any more precise argument to this last claim (which 
is still a conjecture even for the LGW theories). Then a simple calculation consists in 
estimating the "optimal" order n opt beyond which the e-expansion starts to diverge. If 
we only know the first n terms of the expansion, we expect that for n < n opt "ordinary" 
resummation procedures (like Pade) will be sufficient. If n > n opt , or if one seeks higher 
precision, knowledge of the large orders and more sophisticated resummation methods are 
required. Assuming that for e = the n-th coefficient of (3(g) is of order (—Cg) n n\, and 
that we can approximate the fixed point g* by its first order estimate e/B\, the term of 
order n in the e expansion should behave as 



e C 
~B~i 



(70) 



The optimal order n opt is obtained when the absolute value of (70) is the smallest, that 
is for 

n op t 6 - "7J ' ( 71 ) 

where Bi is the one-loop coefficient of the /5-function, and C the large order constant as 
obtained from the instanton calculus. 

With our choice of normalizations for the coupling constant b in the Hamiltonian (4), 
the 1-loop coefficient of the /3-function is 



B 1 = - 
1 2 



(2 - D)S D 
An 



d/2 



r-2 



1 + 



i (r(^)V 
2-d r(^) 



(72) 



with So the volume of the unit sphere in D dimensions 

on = 



r(f) ■ 



(73) 



Let us replace C in Eq. (71) by the variational approximant C var given by Eq. (64). Setting 
finally e = in B\jC (since we are interested in the expansion around e = 0), we obtain 
the following variational estimate for the r.h.s. of (71) 



™ op t e ~ 



16 



(2-D) 



Em 

r(f) 



4 

2-D 



1 + 



i (r(^)y 

2-d r(^) 



(74) 



Let us recall that in practice the e-expansion is used as follows: in order to compute 
for instance the scaling exponent v for a membrane with internal dimension D = 2 in 
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Figure 2: Optimal order n op t(D') for the e-expansion for membrane as function of the extrapolation 
(dimension) parameter D' , as obtained from the variational estimate for the large orders. 

rf-dimensional space, one starts from some point D' ^ D, e = (i.e. d! = AD'/ (2 — D')), 
uses an expansion in e and D — D' (or some more general expansion parameters) to 
evaluate v(D'), which thus depends on the expansion point D'. v is then taken as the 
best estimate u(D' opt ), as determined for instance by a minimal sensitivity criterium. 
Membranes (D = 2) always correspond to e = 4, so setting e = 4 and replacing D by D' 
in Eq. (74), should give an estimate of the "optimal order" n opt (D' ) for the e-expansion 
at D'. The result for n opt (D') is plotted on Fig. 2. 

Some interesting comments can be made on this curve. For D' > 1.6, n opt (D') > 2 
and becomes large as D' — > 2, while for D' < 1.6, n opt (D') < 2 and becomes small as 
D' — > 0. In the first regime (D' — > 2) we thus expect that the power series in e will behave 
like a convergent series, up to some quite large order n opt . In the second regime (D' 
small), we expect that the power series in e will be divergent from the very first terms. 
This is in agreement with the calculations at second order in [20,21]. For large d, the 
2-loop results for v can neatly be resummed, and the stability of the various resummation 
procedures and extrapolation schemes analyzed in [20,21] is good. The final estimates are 
close to the prediction of a variational approximation A/d for v. For smaller values of d 
stability is less good, but in all cases, the reliable extrapolations are obtained for values 
of the extrapolation dimension D' ~ 1.6 or larger. It is not possible to resum safely the 
2-loop results if one starts the e-expansions from D' < 1.5. Thus it seems that our rough 
estimates for the large order behavior may explain some general features of the calculation 
at second order, and corroborate the results of the estimates of [20,21]. 



Of course these arguments are valid if the variational approximation for the instanton 
action stays (at least qualitatively) correct in the limit D — > 2. First let us note that, 
although Eqs. (63) and (64) give estimates for iSi nst and C which are singular when D — > 2, 
our variational formula for n opt is much less singular, since according to Eq. (74) it behaves 



5.3 Limit D 



2 



16 



as 



n opt (D) 



1 16 e" 47 



as D — > 2 , e fixed , 



(75) 



e (2 -/J) 2 

with 7 = 0.577216 the Euler's constant. It has also been noted in [19] that instead of 
using the simple coupling constant b in the Hamiltonian (4) it might be more physical 
and convenient to use as coupling constant the "second virial coefficient" z, defined as 



(2-D)S D 
An 



d/2 



bU 



(76) 



Using z as expansion parameter instead of g = bL € , the large order constant C in Eqs 
(10) and (30) is now 

[2-D)S r 1 ~ d/2 



C 



>D 



An 



(77) 



which in the variational approximation reads 



/7var 



(2 - D)S D 
An 



Cst 



2-D 



as D -> 2 , e fixed . (78) 



(with Cst = |(7re 27 ) 2 £ / 2 . Therefore in this normalization also the singularities as D — > 2 
are simply algebraic. The same remark holds for the coupling constant normalization used 
in [20,21]. 

This prompts us to study the consistency of the variational approximation in the limit 
D — > 2 (or equivalently — > oo) for fixed e. As we shall see, in fact the variational 
approximation becomes exact in that limit. This makes the arguments of this section 
fully valid. 



6 Beyond the variational approximation and 1/d corrections 

In this section we show that the variational result is nothing but the leading term of a 
systematic expansion in 1/d. In ordinary theories like the one described by the effective 
Hamiltonian Tiy which appears in Eq. (50) for a fixed potential V(r) 



H v [r\ = J d D x 



l(v*0 2 + v(f) 



(79) 



this result is not unexpected. For fixed D, the limit d — > oo is nothing but the limit where 
the number of components n = d of the field r becomes large. In this limit, and provided 
that the model is O(n) invariant, it is known that the variational approximation becomes 
exact and a systematic 1/n expansion can be constructed. For the problem considered 
here, the existence of a 1/d expansion is not that easy to prove for two (related) reasons. 
First, the potential V(r) is not fixed, it is also a variable which has to be determined 
self-consistently, and it has a singular behavior when d — > oo. In particular one cannot 
simply take the limit d — > oo while D < 2 is fixed, as can be seen on the exactly solvable 
case of polymers (D = 1) discussed in section 4. Indeed in this case, the equation for the 
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instanton has no physical solution for d > 4. Second, the physically meaningful limit is 
to take d — > oo while e = 2D — d(2 — D)/2 is kept fixed. In this limit, D — > 2 and one 
expects potentially dangerous additional singularities, since the massive free propagator 
is known to have a logarithmic singularity at short distance for D = 2. There is a subtle 
interplay between the corresponding 1/(2 — D) poles and the terms which would naively 
disappear in the large d limit. 



6.1 Expansion around the variational solution 

To evaluate the corrections to the variational approximation, we expand around the vari- 
ational Hamiltonian 



1 2 

2 y 1 2 



(80) 



by writing 

H V [A = Wvarfl - J d D x A(f) (81) 
2 

A(f) = ^£f 2 - \/(f) . (82) 
2 

The saddle point equation (27) which defines the instanton potential V is 

V(f ) + (5(f - r> ))) y = , (83) 
where (• • -)v denotes the expectation value of • • • taken with respect to the Hamiltonian 

for an infinite flat membrane. (Recall that this was the limit we had to take in Eq. (22)). 
This implies that the point Xo can be chosen arbitrarily on the membrane. It is simpler 
to use the Fourier transform of V 

V(k) = J d d fe^V(r) , (85) 

so that Eq. (83) reads 

V-(jfeo) + (e^ xo) ) v = . (86) 

We now expand around the variational solution, using Eq. (81), to rewrite the expectation 
value on the l.h.s. of Eq. (86) as a connected correlation function computed with the 
variational Hamiltonian 

(e^>o)^ = ( e ifcor(*o) . e J* A )^ ar , (87) 

where (• • -) var denotes the expectation value of • • • taken with respect to the variational 
Hamiltonian 7i var 

/ V[r\ . . . e~ nv 

7W\ 



Lvar 



(• ' -)var = J c ^^'-H " ( 88 ) 

/T-var 
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The suffix (■ ■ -) c means the connected correlation function in the usual sense: since 7i var 
is a free Gaussian Hamiltonian, using Wick's theorem, correlation functions like that in 
Eq. (87) can be expressed as Feynman diagrams involving the free variational propagator 
in D dimensions 

d D g p -iq-x 



= / 7^ -TT—T ( 89 ) 



var 



and vertices obtained by expanding the "perturbation term" A(f) of Eq. (82) in powers 
of r; the connected correlation function in Eq. (87) is just given by the restriction to 
connected diagrams. 

Similarly, the free energy density S[V] defined by Eq. (21) can be written as a sum 
over connected vacuum diagrams. 

6.2 Resummation of tadpoles and reorganization in terms of normal products 

In these calculations, we encounter numerous "tadpole diagrams" , which result from the 
evaluation of (r(x)r(y)) var at coinciding points x = y. A standard way to resum these 
tadpoles is to use normal products. This procedure consists in replacing any monomial 
P[r(x)] of the field r(x) at a single point x by the corresponding normal product (or normal 
ordered operator) :V[r(x)} :, defined such that the expectation values of any products of 
normal ordered operators (: V\ : • • • : Vq :) var at non-coinciding points X\ ^ ■ • • ^ Xq 
is equal to the sum over all Feynman diagrams without tadpoles which appear in the 
evaluation of (Vi • • • Vq)^. All normal ordered operators can be obtained 3 from the 
normal ordered exponential, satisfying the relation 

C 2 

where C var is the tadpole diagram amplitude 
with 

C = (47r)- D / 2 r(l- J D/2) . (92) 

C var coincides with the factor of A in Eq. (59) of Sect. 4. Since C var depends on the mass 
m var in Eq. (91) which is chosen to be the mass appearing in 7i V ar, the normal ordered 
products :[•••]: depend explicitly on a mass scale m, and should be denoted by : [• • •] : m . 
This mass scale dependence will be omitted in this section, since we shall always choose 
m = m var . 

With these notations, we can rewrite the operators in Hy and 7i var in terms of normal 
products. This gives 

(r») 2 = dC var l+ :{f{x)f: (93) 
3 by differentiating with respect to k 
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with the identity operator 1. For V(r) we have 

/c\ d k ~ -> 
W? v{k)e ~ 

= £ / S e " ¥cv "^« :( * rT ' : ■ (94) 



If we make the additional assumption that V(r) is rotational invariant, then V(k) depends 
only on \k\ = k and only the even terms m = 2n are non zero in Eq. (94). Integrating 
over k and after some algebra we obtain 



nr~> = £ (-;) 

n=0 V 7 



iV r W 2) At. ^ 



T(n + ci/2) n! 



with the moments Ai n given by 

d d k 



M n = J j^-AkTe-^v(k) 

= 2 ^t£t r dkkd+2n ~ le ~^ C ™V( k ) ■ ( 96 ) 



r(d/2) 

Thus we finally can write Eq. (86) for V(k) as 

V(k) + e-^Cva^.giM.o). e ^ A )^ = 0, (97) 
with the "perturbation" A(r) written in terms of normal products as 

A(f) 



1 + 



1 2 1 A, 

2™™ +2^ 



: r 



?2 . 



n=2 



v( -lV Y{d/2) M n 2 

2. I — ) r(n + d/2) - (r } • • (98) 



Let us first evaluate A(r) when we simply take for V the variational estimate V Yi 
given by Eqs. (58) and (59) 



k 



2 

Cv 



V(k) -> Kar(A;) = -e"^ 1 -™. (99) 
We get for the moments 

M n - Mr = - r( r ( y 2 ) 2) (4vrc var )^ 2 c- (100) 

and we can use the self-consistent equation (60) for m var which is in our notation 

2mLC„ = (47rC var )" d/2 (101) 
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to get for A(f) 

A(r) - A var (f) 



m var C v 



■n—0 \ 



-1 

4C 



2 



: (f») s 



(102) 



The coefficient of the :r 2 : term is zero, since M.\ 



-dm?, 



This means that when we take V W3X as potential V, A is indeed an interaction term, 
which contains no "mass renormalization" , but only higher order interaction terms. 



6.3 A convenient rescaling 

To check whether these terms are unimportant in the large d limit, it is better to rescale 
all quantities in terms of the variational mass m var , that we take as unit of scale. Thus 
let us rescale 

x -> {m vw )~ l x q -> m vai q 

r - {m vw ) D l 2 - l f k - (m vax y- D / 2 k (103) 
V -> (m var ) D V 

In these units, all previous results obtained in this section are given by the same equa- 
tions 4 , provided that we replace 

m var -> 1 , C var -> C , (104) 

that we do not rescale V in Eq. (97), and that we rescale the free energy density £[V] as 

£[V] - m^ r £[V] . (105) 

Thus we got rid of the complicated d and D dependence of m var and keep only the simple 
factor of C given by Eq. (92) in the calculations. 



6.4 The variational solution asD^2 

In order to estimate "how close" the variational potential V yax is from the exact instanton 
potential, let us consider the instanton equation (97), where we replace V by V va _ r and A 
by A var . Then of course Eq. (97) is not satisfied, since 

( :e ifcr(*o). e ^ Avar )^ ar ^ 1 . (106) 

This would be true only if all the n > 2 terms in the expansion (102) for A var were zero. 

We can compute the l.h.s. of Eq. (106) in perturbation theory. The first interaction 
vertices with their coefficients are depicted on Fig. 3. Let us note that the 2-point ver- 
tex, corresponding to the : f 2 : coefficient in Eq. (102), is zero. The "0-point" vertex, 
corresponding to the coefficient of the identity operator 1, disappears in the connected 
correlation functions (since by definition ([•••] l) Cvar = 0), but will be present in the 
free energy density S. The propagator is simply (q 2 + l) -1 . The 0-th order term in the 
expansion of Eq. (106) is 

(:e i£ * xo) :) = 1 (107) 

4 except of course Eq. (101). 
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Figure 3: The interaction vertices from A(f) in the variational approximation. 



because all tadpoles are subtracted by the normal product prescription. On Fig. 4 we 
have depicted the only diagram which appears at order C _1 with its combinatorial weight 
resulting from its symmetry factor, the interaction vertex coefficients and the contractions 
of the d- dimensional indices of the r's. We have also depicted the other possible diagrams 
which do not contribute, since they contain a tadpole which is subtracted by the normal 
order prescription. The diagrams at order C~ 2 are depicted on Fig. 5, together with their 
weight. In general, these diagrams contain internal closed loops, which give a factor of d, 
and open chains which must end at some kr(x ) in the exponential :exp (ikf(x )) :, thus 
giving a factor of k 2 . 

We can now look at the limit when the internal dimension of the membrane D goes 
to 2, while e is fixed. In this limit, the bulk dimension of space d goes to infinity as 



d = 



2(2D - e) 
2-D 



8 -2e 
2-D 



:io8i 



and the tadpole amplitude C given by Eq. (92) diverges like d, since 



C 



(Au)- D I 2 Y{\- D/2) ~ 



1 



2ir 2-D 



~ d 



1 



4tt(4 - e) 



(109) 



On the other hand, the Feynman amplitude of the first diagram depicted on Fig. 4 and 
of those of Fig. 5 are finite when D — > 2, since they do not have any long distance (infra- 
red) or short distance (ultra-violet) divergences for < D < 2. This implies that as 
D — > 2, the contributions of diagrams that we are considering vanish at least as fast as 
2 - D ~ l/d. 

This result is in fact valid for all the diagrams which appear at higher orders in the 
evaluation of the l.h.s. of Eq. (106). All diagrams vanish individually at least as l/d, 
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Figure 4: First order graphs in the variational approximation. 
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Figure 5: Second order graphs in the variational approximation. 

and it is possible to generate a systematic expansion in powers of 1/d. Of course there 
will be an infinite number of diagrams which contribute at a given order in 1/d, that we 
shall characterize later. Let us assume that the sum of all diagrams which contribute 
at a given order is convergent. (This turns out to be true at least as long as e > 0.) 
Then this implies that the saddle point equation (97) holds at leading order in 1/d for 
the variational solution 

V vai {k) + e"^ c ™ (:e^ (xo) : e^ A ™)var = e~^ c ™ , (110) 

when we take the limit (2 — D) ~ d^ 1 — > 0, e and k fixed. 

A similar statement can be made for the vacuum diagrams which appear in the per- 
turbative expansion of the free energy density. The first terms of this expansion are 
depicted on Fig. 6. The first two graphs denote the tree and one-loop contributions given 
respectively by (minus) the first term in the expansion of A var 

Sffi = " C ( 2 + ^) (I") 

and by the loop integral 

<8 = !/§>«Nfl^c. (112) 
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Figure 6: Free energy density. 
The sum of these two terms gives the variational free energy density 

£var = C (-2~ + £) = (113) 

which is of order 0(d), while higher order corrections are at least of order 0(1) (in our 
normalization where m var = 1). 

What are the consequences of this remarkable fact? First it will be possible to solve 
the saddle point equation (97) order by order in 1/d, and to show that the exact instanton 
potential V differs from the variational solution only by corrections of order 1/d, 

v(k) = K ar (*0 [1 + 0(01 ( 114 ) 

which are finite as long as e > 0. Using this fact we can show that the variational instanton 
action, which is in our normalizations 

Svar = £var + \ J ^ = (} ~ J}) C > ( U5 ) 

is of order 0(d), and differs from the exact instanton action by subdominant terms of 
order 0(1). This justifies the use of the variational method as well as the large order 
analysis of the e-expansion results made in the previous Section. 



6.5 Leading 1/d correction for the instanton potential 

In order to study the 1/d corrections, let us rewrite the normal product expansion (95) 
for the exact instanton potential V as 



n=0 



With our rescaling in Eq. (103), starting from Eq. (96), the moments /J, n are given by 

o rnn+d/2 poo 

= v( \, l0 , / dkk^e-^V(k)- (H7) 
r(ra + d/2) J 

We assume that V differs from V var by 0(d~ l ), or equivalently that 

Hn = "I + 5 ~J , $n = 0(1) . (118) 
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Figure 7: The vertices from A(f), Eq. (119). 

The perturbation A(r) given by Eq. (98) is then 

m = C g - 2„) + * :f : + (-2C) g g)" £ . (119) 

This generates the 2n-point interaction vertices for the perturbative expansion around 
the variational Hamiltonian depicted on Fig. 7. 

The first vertex with n = is just the correction to the free energy. The additional 
2-point vertex is a mass correction and is of order \jd. The 2n-point vertices (n > 2) are 
similar to those of Fig. 3 and are of order 1/d™ -1 . 

Before embarking onto the detailed calculations, let us recall what we are about to 
do. Our equation for the exact instanton potential V is a self-consistent equation (since 
both sides contain V), and was written as a perturbative expansion with respect to the 
variational solution as 

-V(k)e Ck2/2 = (:e ikr ^:eI* A ) C . (120) 

\ / var 

A, given in Eq. (119), contains all terms of the exact potential V with the exception of 
those already taken into account in the variational Hamiltonian. 

It is then easy to see that the perturbative expansion generates diagrams, which can be 
organized in a 1/d expansion in terms of the chain of bubbles depicted on Fig. 8. Indeed, 
each 4-point vertex carries a factor of C ~ 1/d and each bubble carries a factor of d, so 
that the whole chain is of order 1/d. A careful but simple analysis shows that only four 
different classes of diagrams contribute to the r.h.s. of Eq. (120) at order 1/d; they are 
depicted on Fig. 9. The suffix Yn \ (n = 0, 1) refers to the minimal number of bubbles in 



the chain in order not to have tadpoles. The situation is thus quite similar to the 1/n 
expansion in models with a n-component field with 0(n) symmetry. 
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Figure 8: The chain of bubbles 
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Figure 9: 1/d diagrams for the expansion of (:e lkr(x °' ) : e-^ Avar )^, 



The amplitude for a single bubble diagram with external momentum p is 

d D q 1 1 r(2-D/2) 



) D g 2 + 1 (p + g) 2 + 1 



( 47r )D/2 



■J(P) 



:i2r 



with 



J(p) = / drr [l + :r(l -a;)p 2 p 2 (122) 
Jo 

and can be expressed in terms of hyp ergeo metric functions. In particular one has 

= (1-D/2)C . (123) 



Taking into account the symmetry factor, the amplitude for the [n] -truncated chain of 
Fig. 8 is the geometric function 



n 



= HW(p) = £ 



m>n 



(124) 



so that in particular the untruncated chain is 







M(°)(p) 



1 + A*2 



2D - e 



J(p) 



(125) 



Similarly, the free energy density £ can be expanded in 1/d in terms of diagrams 
involving chains of bubbles. The diagrams which are present at order 0(1) are depicted 
on Fig. 10 



The last diagram is the [2j -truncated closed chain of bubbles, whose amplitude is 
(taking into account the symmetry factors) 



I 
I 

* 



n 



', = G<"> = 



E 1 

< J m 



d D p 

(2ir) D in 

m>n 



-^)^B(p) 



(126) 
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Figure 10: 1/d diagrams for the expansion of the free energy density £ at order 0(1) 

With these notations, we now compute the 1/d correction for V. With the symmetry 
factors for the diagrams of Fig. 9, Eq. (97) for V is at that order 




(127) 



Since we are interested in computing V at order 1/d only, and since the diagrams on the 
r.h.s. of Eq. (127) are of order 0{1), we can use the ansatz (118) and replace /i 2 and /X3 
by —1 on the r.h.s. of (127) (but keep 5\). Now we use Eq. (123), which implies that 




(128) 



and from Eq. (124) we have 




(129) 
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We can also rewrite 



We obtain 



V(k) 



+ 




4C 



t \-d/J,2 



% * 



-fc 2 C/2 



1 + (-A; 2 ) 
\ 



2(4-e)C 

+ 



V 



e-4 
2de 



C 



(130) 




+ (^ 4 ) : a ; -» + o 



;i3i) 



Now we have to compute Si. It is fixed by Eq. (117) which relates the /i n 's to V. Indeed 
we have from Eqs. (117) and (118) 

9 (pl+d/2 roc 

S ^« + d T(^W)l (132) 

and using Eq. (131), performing the fc-integral and keeping only the leading terms as 
d — > oo, we obtain the equation 




Finally, we can use the simple relation (derived in appendix B) 



mi = (4-^)4 



i + 




(134) 



to express Si in terms of the simple closed chain of bubbles defined by Eq. (126) 

+ o(K\ • (135) 



Si = 



d 



2C(2-e) V 
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Inserting this result into Eq. (131), we obtain the potential V, and we can calculate the 
moments //„ from Eq. (117). Since for n <C d, the moments are independent of n, we find 
that 

^ = -1 + 5 j- + O (J-^J Vn > . (136) 

Therefore the ansatz (118) is consistent and our calculation at order 1/d makes sense, 
provided that the closed chain of bubbles G^ 2 \ which is a function of D, e and /i 2 , has a 
finite limit as D — > 2 (and /i 2 = —1). This is the case as long as e > 0. Indeed, the chain 
G (2) is given by Eq. (126) and using Eqs. (121) and (122), it is given in the limit D — > 2 
by 



G{e) = : 



D=2 



d 2 q 



- log 



1 - 



i-i)jfo)] ■ (l-j)'(v) 



(137) 



where 



■/(?) 



HQ)\D=2 

2 



/ dx [l + x(l — x)q 2 ^ 1 



9 vV + 4 



log 



yV + 4 + g 
vV + 4 - g 



138) 



Indeed, we have for J(q) 



< J(q) < 1 Wq ^ ; J(g) ~ 1 - as q -> ; J(g) ~ — log(g 2 ) as g -> oo , (139) 



hence the integral (137) defining the function G(e) is convergent for any e > 0. 
Similarly, the chain with a marked point is given in the limit D — > 2 by 



G(e) 



D=2 



/ 



d 2 g 

(27T) 2 



(!-!)■/(?) 1 ^ 4 



■/(«) 



- (4-e)-^). 



(140) 



It is finite for e > 0, but has a logarithmic divergence at e = 0, since from (139) the 
integrant in (140) behaves as g~ 2 as q — > when e = 0. A simple calculation shows that 
G has a logarithmic singularity of the form 



G{e) ~ |^log n 



as 







(141) 
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With these notations, the final result for the instanton potential is at order 1/d 



V(k) 



,-fc 2 C/2 



1 



d V V2-e 



-GW + i^G'( £ ))-^( £ )) + o(i 



(142) 



We note that the corrections to the variational potential are finite as long as < e < 2, 
but are singular at e = and e = 2. The singularity at e = 2 is not surprising. We 
have seen from Eq. (61) that the variational mass m var itself is singular at e = D, which 
corresponds in the limit D — > 2 to e = 2. Therefore it is expected that the corrections to 
the variational result will also be singular. The singularity at e = is new, since m var is 
regular at e = 0. We shall come back to this issue later. 



6.6 1/d correction to the instanton action 

It is now easy to compute the corrections to the instanton action S, defined by Eq. (24). 
The free energy density S has the following graphical expansion 



S(V) = -C (^-2^ 



+ 



d 




■ 



+ O 



(143) 



where the amplitude for the simple loop is given by Eq. (112) 




/ 



d D q 



log [1 + 1/q 2 ] 



D 



C 



(144) 



The second term in Eq. (24) is 

1 



1 



2 , dV V(f? = 2 



d d k 
(2tt)° 



\V(k)f 



Within our rescaling (103) this term is 



5 I ! " 



c 



2C d/2 



dk k d - l \V{k)\< 



r(d/2) Jo 

and from Eqs. (117) and (118) it is simply obtained from the moment /i by 



v 2 



c 



-l-2^o +0{d- 2 ) 



(145) 



(146) 



(147) 



Thus the instanton action, which is the sum of Eqs. (143) and (147) is finally given at 
order 1/d by 



S(V) = C [ d - + A 




1/ 
2 1 



+ ° { d 



(148) 
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= C(l--i) - 1 -G^(D,e)+o(^j . (149) 

This result holds within the rescaling (103). The first term on the r.h.s. of Eq. (149) is 
the variational result, and is of order 0(d). The second term is of order 0(1). Thus the 
first correction to the variational result is given by the closed loop diagram of Eq. (137). 

According to Eq. (105), to recover the normalization used in the previous sections, we 
have simply to multiply Eq. (149) by m® ar , where m var is the variational mass as given by 
Eq. (61). 

6.7 Discussion of the result 

(1) The corrections of the variational instanton potential do in fact not contribute to the 
correction of the instanton action. Indeed, they are entirely contained in the moment /j,q 
which disappears in Eq. (149). This is not surprising, since the instanton potential V is 
defined by a variational principle with respect to V. 

(2) The correction to the instanton action is finite as e — > 0, while from Eq. (142) the 
correction to the potential V diverges as log(e). The divergences occur in fact because 
the chain of bubbles of Fig. 8, M^(q), has a singularity at small momentum q as D — > 2 
and e — > 0. According to Eqs. (125) and (139), it behaves as 

1 

" "© " ' ~ q2/6 + e/4 + 0{2 ~ D) ' 

In other word, this chain of bubbles behaves as a massless propagator at e = and D = 2. 
We have no physical interpretation of these IR singularities, and why they cancel in the 
instanton action. However, the fact that as D — > 2 the instanton potential is singular as 
e — > is similar to the situation at D — 1 discussed in section 5.1. Here also, we know 
from the exact solution that the instanton potential is singular when e — > (d — > 4), 
while the instanton action is regular, and in fact Eq. (67) corresponds also to a logarith- 
mic behavior as e — > 0. Therefore, it is reasonable to conjecture that at next orders in 
1/d, the IR singularities which occur as e — > in the instanton potential still disappear 
in the instanton action. This would imply that the coefficient C has a regular behavior 
when e — > and that the limits d — > oo and e — > can be exchanged, providing further 
consistency to our arguments for the large order behavior of the e-expansions for the SAM 
model. 

(3) The 1/d correction to the instanton action iS inst is negative, as expected, since it 
should improve the variational estimate for <Si nst , which is an upper bound. 

(4) It is interesting to see if the 1/d correction can be used to improve the variational 
estimates for quantities of physical significance, such as the large order constant C. Let 
us start from the result (149) for the instanton action S at first order in 1/d, and insert 
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Figure 11: The function G(e) denned in Eq. (137) from the chain of bubbles in D = 2. 



it into Eq. (30) for C. We obtain 

^ War 



2C{D) 



(151) 



Now we only keep the leading l/d corrections on the r.h.s. of Eq. (151) in the limit d — > oo, 
for fixed e by using 



G^(D,e) = G(e) + 0(O , C(£>) = A J_ + 0(1) , 

4 7T 4 — e 

with G(e) given by (137). We thus obtain 



C _1 = C' 1 

u War 



1 - 2tt(4 - e)G(e) - + C(rf~ 2 ) 
o 



(152) 



(153) 



We have plotted on Fig. 11 the function G(e) for < e < 4, as obtained by a straightfor- 
ward numerical integration. G is maximal for e = 0, where we have G(0) = 0.6014 

Let us estimate the first l/d correction for the case of polymers in d = 4 dimensions. 
We simply have to set d = 4 and e = in (153). We find that the l/d correction is 



2vr(4- e )G( e )- 



2ttG(0) = 3.78. 



(154) 



which is much larger than 1 ! Thus the l/d correction is very large for d — 4 and one 
should take into account the subleading corrections and resum them. We recall that for 

d = 4 and e = we have C~^ act /C~^. = 2/3 = 0.666 In practice, we expect that 

for e = the l/d corrections will be smaller or equal to the variational O(d ) result for 
d > 16 ! Therefore we cannot use naively the l/d corrections that we have calculated 
to improve the variational estimates. In the next subsection we propose an improved 
resummation procedure, which takes into account some of the higher order corrections 
which are contained in the bubble chain diagrams, and which gives much better results. 
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6.8 An attempt to go beyond the first 1/d correction 

In the last section we have shown that a straightforward 1/d-expansion cannot be applied 
to small dimension, as e.g. polymers in d = 4. In this section, we propose a different 
approximation scheme. It consists in summing exactly the chain of bubbles G^ 2 -*, and 
keeping the full D, (/-dependence when extrapolating to low dimensions, instead of ex- 
panding this quantity about d = oo. First, in Eq. (126), only the leading (/-dependence 
had been kept. Taking into account the complete (/-dependence, we obtain 



* 

* 
■ 



n 



(n) 



d + 2 



E 1 

< J rn 



dp 

(2n) D — //; 

m>n 



-^2) 



4C 



(155) 



Eliminating as before M(p) in favor of JJ(p), d in favor of e and replacing /i 2 by its leading 
contribution fj, 2 = —1 (justified later), we obtain for the correction to the variational 
result 



S-S v 



D{2D-e) sin^ 
(2 + D-e)(D-e) n 



x 



dpp*- 1 ln^l 



2 + A 2 + D-e,, ' 

— - A — J(p) + — - A — Hp) 



(156) 



In the remainder, we shall focus on the case D — 1, for which we can most easily test Eq. 
(156). In D = 1, J(p) is exactly given by 



4+p 2 ■ 

Eq. (156) is then integrated (using the residue calculus) with the result 

2d 



S-S v 



D=l 




(d + 2)(d-2) 



The large order estimate is finally obtained as 

d 



C~\D = 1) ~ 2n d/2 



1 + 



d + 2 



d d 7 
1 "2-8-4 1 



(157) 



(158) 



(159) 



which is plotted on Fig. 12. We see that this corrects 50% of the deviation of the variational 
result from the exact result in d = 4, and is even better in lower dimensions. Note that 
this is not the straightforward l/rf-correction, obtained in section 6.7. In the following, 
we want to justify that the result given in Eqs. (158) and (159) is meaningful for small 
dimension, by analyzing the case D = 1 and d — 4, i.e. e = 0, where calculations are most 
easily done. 

The first observation is that rf ) (p) is nicely convergent for all values of p. In contrast 
to D = 2, no divergence at p = appears. 
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Figure 12: The inverse of the large order constant l/C for the Edwards model (D = 1) as a function of 
the bulk dimension d. The dotted curve is the 0(1) variational estimate, the dashed curve the estimate 
from Eq. (159), the continuous curve the exact result. 

The next observation is that all terms in the chain of bubbles G^ 2 -* are positive and the 
sum rapidly converging. (If we denote by 100% the difference between variational estimate 
and exact result, the contributions of the first terms are 28%, 11%, 4%,. . . , respectively, 
which add up to the 50% given above.) 

We also observe that most of the missing contributions to S(V) are positive: Let us 
call a vertex even, if it possesses an even number of pairs of lines, and odd otherwise. Any 
even vertex then contributes a factor of 1, whereas any odd vertex contributes a factor 
of (—1) (see Eqs. (116) and (118)). Since all integrals are positive, in order to build up 
a negative contribution, an odd number of odd vertices has to be taken. One sees by 
inspection, that this is more difficult than for a positive diagram (containing even vertices 
and an even number of odd vertices). The first diagram of this kind is 




(160) 



Another hint comes from a second class of diagrams, which can be summed, namely 
all "watermelon" diagrams with 2 vertices only: 




(161) 
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The first diagram, equivalent to 28%, is already contained in <&( 2 \ the following contribute 
6%, 1.5%, . . . respectively. They can be resummed and are then given by the integral 




(162) 



which contributes about 8.2% in d — 4, when neglecting the first diagram already taken 
into account in G^. 

A more systematic approach to understand why higher order vertices are subdominant, 
is to look for an additional small parameter. We have seen that C ~ d for d — > oo. 
Comparing the ratio d/C at d = oo and d = 4, we obtain (always for e = 0) 

d/C\ d=A _ 1 
d/C\ d=00 2tt 

This suggests (as exemplified by the preceding calculation) that higher order vertices, 
which come with additional factors of 1/C, indeed contribute an additional small factor. 

The last point to verify, is that also the l/rf-correction to the potential is small. In 
the notation of Eqs. (116) and (118), the parameter Si/d, given in Eq. (135), reads 




* - * £>=l,e=0 

It is indeed small. 

Another important limit to verify is d — > 0. We have argued in section 5.1 that in 
that limit the variational result becomes exact. Our correction in Eq. (156) consistently 
vanishes. Diagramatically this can be seen from the fact that any diagram which con- 
tributes to the free energy has at least one closed loop, and therefore comes with a factor 
of d. This is also the leading contribution for small d: due to the "islands and bridges 
theorem", which states that if from any island an even number of bridges starts, it is 
always possible to construct a path which uses any bridge exactly once, the minimal 
number of closed loops, and thus factors of d, for any given diagram is one. It is therefore 
impossible to simply reorganize the perturbative expansion in a <i-expansion. However, 
the large-<i expansion should respect this exact property. This is not assured by the naive 
l/o?-expansion, but by our modified result discussed here. 

All these arguments suggest that Eq. (159) is a sensible correction to the variational 
result, and should be improved by taking into account higher order vertices, although this 
expansion is not a systematic expansion in a small parameter. 

7 Conclusion 

Let us first summarize the results of this paper. We have developed a general formulation 
to estimate the large orders of perturbation theory for the Edwards model of self-avoiding 
membranes and polymers. We have shown that these large orders are controlled by 

35 



(163) 



a semi-classical effective potential V(r), solution of a non-local extremization problem. 
This effective potential is the analog to the Lipatov instanton in Landau- Ginsburg- Wilson 
models. In the case of polymers (membranes with internal dimension D — 1), this SAM- 
instanton corresponds precisely to the Lipatov instanton for the n = components LGW 
model. The large order behavior for the SAM model is derived in the general case (D ^ 1). 
The equations for the SAM-instanton are solved within a Gaussian variational approxi- 
mation. The result is for D — 1 compared to the exact results from Lipatov's method, 
and is found to be qualitatively correct for polymers in d = 4 dimensions. Finally it is 
shown that the variational result is the first term of a systematic expansion in 1/d, where 
d is the dimension of space in which the membrane fluctuates. 

All these results are new, and represent a considerable advance in the understanding 
of self-avoiding membranes beyond the first orders of perturbation theory. A number of 
interesting issues still has to be addressed: 

- We have seen that a systematic expansion in 1/d can be constructed, of which the vari- 
ational approximation is only the leading term. However, adding the first correction in 
1/d to the instanton action does not give reliable results for small d. This is a numerical 
problem which comes from the fact that in the interesting cases, the corrections in 1/d are 
already quite large. We have proposed a modified procedure to resum the first corrections 
which gives better results when applied to the case of polymers (D — 1, d — 4), but more 
systematic resummation methods are needed in order to improve in a reliable way the 
variational results. 

- We have not calculated the contributions from the fluctuations around the instanton, 
which should give the value of the global constant A' in Eq. (32). This calculation is 
technically more difficult than the one for the instanton action itself. 

- We have only discussed shortly and at a qualitative level the consequences of our re- 
sults for the e-expansion of the scaling exponents. It would be interesting to obtain more 
precise results. In particular, the discussion relies on the limit e — > for the large order 
estimates. We have shown that this limit exists in the variational approximation, but 
that the 1/d corrections then suffer from IR divergences. These divergences cancel for 
the instanton action at first order in 1/d, and we conjectured that this feature persists at 
higher orders. Further studies are needed to clarify this important problem. 

- Finally, let us mention that our approach can be applied to other classes of interactions 
between membranes. 
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A The <£> 4 instanton in the limit d — > 4 

The instanton equation (37) for \l/ is, assuming rotational invariance i.e. ^o(f) = \& ( r ) 
with r — \r\, 

% + + 2E ^ + 2^1 = . (165) 

We look for the solution which is regular at the origin, which has no zero for finite r and 
which vanishes at infinity, i.e. 

< x I / o( r ) < oo for < r < oo , ^o(r) — > when r — > oo . (166) 

£"0 < is fixed by the normalization (38) 



i*oir 



y dV^rf = — — — drr^^r) 2 = 1. (167) 



From Eq. (45) we know that E can be extracted from the action of the instanton for the 
LGW action (44), which is finite when d — > 4, such that _E behaves as 

£ oc d-A . (168) 

For d = 4 and _E = 0, the general solution of Eq. (165) is not normalizable and reads 

*o(r) = ^ • (169) 

For d 7^ 4, E < and r — > 00, \[> is a solution of the linearized equation \P" + (d — 
l)/r^' + 2Eq^> = 0, given by a Bessel function 

* (r) ^ ^-^^/^(^Sol^r) . (170) 

Assuming that 

r 2 I S 1 < 1 when d -> 4 , (171) 

it has the asymptotics 

f r 2 ~ d for r 2 |2£ | < 1 , , 

^o{r) oc < i-d i9p„ii/2_ r , lnn , . (172) 

v 7 [ r — e-l 2 ^! 7 r for r 2 \2E \ > 1 v 7 

For r C r <C |2_E | _1 ^ 2 ! ^0 decays algebraically as 1/r 2 . It is exponentially decreasing 
for r (— 2E )~ 1 / 2 . We may thus approximate the integral in Eq. (167) by 



a 

*o|| 2 ~ 2tt 2 y^ drr 3 « 2tt 2 r 2 log 



W4 - d. 



;i73) 



where a is the proportionality factor in Eq. (168). From the constraint \\^o\\ = 1 we 
deduce that the size of the instanton r becomes small as d — > 4 and scales as 

r » ~ ■ (174) 

This is consistent with the assumption (171). 
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B A simple relation 



We derive the relation 



mi = (4-i>)4 




;i?5) 



For that purpose, we write the amplitude for the last diagram as an integral over p, the 
D-dimensional momentum which flows through the chain of bubbles 



® 



dD P D(p)eW(p) = 2(47r 



(27T) 



r(D/2) 



_D/2 /-co 

/ dpp D_1 D(p)HW(p) . (176) 



The definition (124) for the chain of bubbles was 



H (1) (p) 



i +/i 2 



4C 



- 1 



and B(p) is the amplitude for the bubble with a single mass insertion 

d D q 1 1 



D(p) = p 




(2tt)^ (g 2 + l) 2 (p + g~) 2 + 1 ' 



;i77) 



(178) 



We are careful in distinguishing the vector p from its modulus p = |pj. Similarly the 
amplitudes for the first two diagrams are 
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and 
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2(4tt) 
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l + /i 2 ^B(p) 



1 + /i2 4C 



We can obtain D(p) from B(p). Let us introduce a mass m and consider 



/" d 5 1 



(2ir) D (,f + m 2 ) 2 (p + g*) 2 + m 2 
and similarly let us introduce a mass in the bubble B of Eq. (121) and define 



B(p, m) 



1 



d^g 1 
(27r) D ((f 2 + m 2 ) (p + g*) 2 + m 2 
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(180) 
(181) 

(182) 



It is easy to see that one has 



to- — B(p, to) = -4m 2 D(p,m) , 
dm 

and since by homogeneity 

B(p,m) = m D - 4 M(p/m) , 
we deduce that (setting at the end m = 1) 

(D-4)B(p) - p^-B(p) = -4D(p). 
dp 



(183) 
(184) 

(185) 



We can use Eq. (185) to express D(p) in terms of B(p) and its derivative in Eq. (176) and 
integrate by part in order to eliminate the derivative of B. We obtain 
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dpp 



D-l 



1 - 



D 
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Using Eqs. (179) and (180) we obtain the identity (175). 



(186) 
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